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We examine the real space renormalization group method of finding excited eigenstate (RSRG-X) 
of the XX spin-1/2 chain, from entanglement perspectives. Eigenmodes of entanglement Hamilto¬ 
nian, especially the maximally entangled mode and corresponding entanglement energies are studied 
and compared with predictions of RSRG-X. Our numerical results demonstrate the accuracy of the 
RSRG-X method in the strong disorder limit, and quantify its error when applied to weak disor¬ 
der regime. Overall, our results validate the RSRG-X method qualitatively, but also show that its 
accuracy decreases with increasing temperature. 

PACS numbers: 


I. INTRODUCTION 

Entanglement has emerged as an increasingly impor¬ 
tant way to characterize phases and phase transitions 
in quantum many-particle systems^. Among different 
measures of entanglement, von Neumann entanglement 
entropy (EE) remains the most popular one, and it is 
defined as follows. If a system is in a pure state |'I'), 
the density matrix of the system is p = Iff) ('I'|; EE is 
the von Neumann entropy of the reduced density matrix 
(RDM) of a chosen subsystem, which in turn is obtained 
by tracing over degrees of freedom outside of the subsys¬ 
tem. Since RDM is a positive definite operator, we can 
write the RDM of subsystem A, as exponential of a 
Hermitian operator, H^: 


where Z = Tr{e~^ ), and we call entanglement 
Hamiltonian of subsystem A. Entanglement spectrum 
(ES), which is the spectrum of has also been found 
to be highly valuable^. 

In addition to the EE and ES, we recently studied 
eigenmodes of entanglement Hamiltonian. For the spe¬ 
cial case of free fermion system (that we consider in this 
and earlier papers), the entanglement Hamiltonian 
is also a free fermion Hamiltoniani^ 

Na 

( 2 ) 

bi=i 

We find the single-particle eigenmode that contributes 
most to EE, the maximally entangled mode (MEM) 
(see below for definition), contains rich pWsical infor¬ 
mation about the system^i^. In Ref. |3, we stud¬ 
ied the ground state entanglement properties of random 
anti-ferromagnetic XX spin chain and compared it with 
the predictions of the real space renormalization group 
(RSRG) method^. At each step of RSRG method, two 
spins form a singlet and the ground state is approximated 


by tensor product of singlets formed with all lengths. Ex¬ 
pectations of RSRG about entanglement are that MEM 
be highly localized at the location of two spins that 
are entangled and also entanglement energies be either 
0 for such entangled spins that cross the boundary, or 
±oo otherwise. Our study showed that these predic¬ 
tions of RSRG on ground state are indeed correct in the 
strong disorder regime. Furthermore, in the weak disor¬ 
der regime, where RSRG remains qualitatively correct at 
long length scales but quantitatively inaccurate at short 
distances, MEM gives the profiles of the effective spins 
that form singlets. 

Recently more attention has been focused on entangle¬ 
ment entropy of highly excited states (namely states with 
finite excitation energy density), including free fermion 
systemsi^— . In the meantime, a new study extended 
the idea of RSRG to include calculation of excited eigen¬ 
state of random spin chains, which is named RSRG-X 
methodii. This method relies on the fact that at each 
step of eliminating spins coupled by the strongest bond, 
one may choose to put them in the ground state or one of 
the excited states of the bond. For example, in the case 
of XX spin chain, either singlet or each one of the triplet 
states can be chosen for two spins, based on probabilities 
given by the temperature. Repeating this procedure gen¬ 
erates a typical highly excited state for a given tempera¬ 
ture. Justification of RSRG-X is subtler than RSRG, as 
there exist distinctive excited states generated by RSRG- 
X whose energies are very close, which may be mixed by 
interactions that are left out in RSRG-X. 

The main purpose of this paper is to critically exam¬ 
ine the RSRG-X method using insights from entangle¬ 
ment, and establish its validity by comparing its pre¬ 
dictions about the entanglement properties of the sys¬ 
tem with exact numerical calculations. We use meth¬ 
ods similar to Ref. 0, with substantial extensions needed 
to handle highly excited states. As detailed in the fol¬ 
lowing sections, our results agree with the predictions 
of RSRG-X method qualitatively. We are also able to 
quantify the error RSRG-X makes, and show that gener¬ 
ally speaking such error increases with increasing energy 
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density/temperature, and decreasing disorder strength. 

This paper is organized as follows. In sec. [iTl the 
model that we study and also terminology about EE 
that is used in this paper and calculation methods are 
introduced. Sec. Ell is devoted to the explaining of the 
RSRG/RSRG-X method and applying it to our model. 
The main section is IIVI which has two subsections, in one 
we study entangled modes and the other is about entan¬ 
glement energy. Einally, a summary is given in the sec. 
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eigenvalues of correlation function, nf, and eigenvalues 
of h^, €^, are related to each other: 


nk = 


I -I- 


( 8 ) 


where referred to as entanglement energy, and is 
the probability of finding a fermion in the corresponding 
mode. In terms of the eigenvalues, the entanglement 
entropy is: 


II. MODEL AND ENTANGLED MODES 

The model we work with is a one-dimensional (ID) 
spin-1/2 XX model with N sites and with random nearest 
neighbor couplings J. The Hamiltonian of the system is 

N-l 

n—1 

We use open boundary condition. By mapping spin op¬ 
erators to fermion operators via Jordan-Wigner transfor¬ 
mation: 

iTrW.. si's” — / A\ 

Cn = e ^ , (4) 

where c(c^) is fermion annihilation (creation) operator, 
the fermionic representation of the Hamiltonian is: 

I w-i 

^ ~ 2 E! Jn{c\Cn+l + c]^_|_iC„). (5) 

n—1 

J’s in Eq. @ will be generated based on random dis¬ 
tribution functions to be specified later. The number of 
fermions is related to magnetization as: 

1 ^ 

+ (6) 

n—1 

In this work we diagonalize Eq. m and generate its 
exact eigenstates, study their entanglement properties, 
and compare with predictions of RSRG-X. 

To study bipartite EE, we divide the system into two 
subsystems A and B, often (but not necessarily) with 
equal number of sites. Subsystem A is from site 1 to Na, 
and subsystem B is from site Na -I- 1 to V. RDM of 
subsystem A is written as Eq. ([T]) which contains entan¬ 
glement Hamiltonian written as Eq. ©• To determine 
eigenmodes and eigenvalues of we follow Ref. [l^ by 
defining correlation function 

=< cjc, >, (7) 

which is an Na x Na Hermitian matrix. Eigenmodes of 
the correlation function are the same as those of h^, while 


Nf 

EE = -'^[n^ In(n^) + (1 - n^) ln(l - n^)]. (9) 

fc=i 

There is no contribution to the EE when = 1 {n^ = 
0) which corresponds to = —oo (e^ = -too) and the 
most contributions comes from n^’s close to 1/2 (i.e. 
be close to 0). 

RDM of subsystem B can also be written as the ex¬ 
ponential of a Hermitian operator similar to Eq. o, 
and also a correlation function for subsystem B can be 
defined similar to Eq. ©• Thus, entanglement Hamil¬ 
tonian eigenmodes of both subsystem A and B can be 
obtained. We use the Klich methodic to stick together 
eigenmodes of subsystem A (which are from site 1 to 
Na) and eigenmodes of subsystem B (which are from site 
Vyi-l-1 to N). Projecting Klich eigenmodes to the subsys¬ 
tem A (B) yields to the entanglement eigenmodes of the 
subsystem A [B) weighted by the probability of finding a 
fermion in those eigenmodes. For each Klich eigenmode, 
which is filled by one fermion, -I- nf = 1 so = —ef. 
Those Klich eigenmodes that correspond to small mag¬ 
nitude of e (and thus contribute the most to the EE) are 
of particular interest. In our previous paper— we showed 
that the Klich eigenmode corresponding to the smallest 
magnitude of e (which we call maximally entangled mode 
(MEM)) contains very useful physical information about 
the system. 

As a simple example, let us consider the case of two 
spins. We choose one of them to be subsystem A and 
the other subsystem B. In the corresponding fermionic 
representation there are two energy levels. Two spins 
can be in each one of four different states: singlet or 
one of three triplet states. If the two spins are in 
triplet^^=|4,4,), then by Eq. ([SI), Np = 0 and thus 
r^A _ Qj. _ _|_QQ^ jf spins 

are in triplet-f-i-=|tt)j where Np = 2 and = 1, 

or = — oo. Obviously, the two spins are not 

entangled in these cases. But when two spins are in 
singlet=-^(|t4') — Ij-t)) (where there is one fermion in the 

lower energy level) or tripletf4,=-^(|ti) + lit)) (where 
there is one fermion in the higher energy level) then 
Np = 1 and = 1/2, and two spins are en¬ 

tangled. 
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TABLE I: Eigenstates and eigenvalues of two spins, the cor¬ 
responding probability associated with temperature T, and 
effective spins in the RSRG-X method. Z = 2 + 2 cosh ^. 


Eigenstate 

Eigenvalue 

Probability 

Effective 

coupling 

singlet^ ^[|t 4 > - lit)] 

-J/2 

1 J/2T 

J « 

triplet4.4=^[|ti) + |4-t)] 

+J/2 

1 ^-J/2T 


triplet 4 . 4 .= |4t) 

0 

1 

z 


triplet 44 = U) 

0 

1 

z 

J 


III. RSRG AND RSRG-X APPLIED TO 
RANDOM XX MODEL 

Consider anti-ferromagnetic XX spin-1/2 chain with 
coupling constants J’s that are randomly distributed by a 
distribution function. The approximate ground state can 
be obtained using RSRG method as following: at each 
step of RSRG we choose the biggest J and consider two 
spins that are coupled by this Jmax (spins number 2 and 

3 in Fig. [T|). We put these two spins to be in the lowest 
energy state, the singlet state. Since coupling of these 
two spins is bigger than neighboring couplings Jl and Jr, 
we treat them as perturbations. First order perturbation 
vanishes and by using the second order perturbation we 
obtain an effective coupling between spins number 1 and 

4 to be J ~ and we remove spins number 2 and 3. 
This procedure continues until we form N/2 singlet pairs 
(assuming even N). The ground state of the system is 
approximated as the direct product of these singlet states 
within RSRG. 

Excited eigenstates can also be obtained by a modified 
version of RSRG which is called RSRG-Xii. In this mod¬ 
ified method, two spins with the largest coupling constant 
are chosen to be in one of the four eigenstates with energy 
E by the Boltzmann distribution: 

P = Z = 2 + 2 cosh(10) 

namely for each one of the singlet/triplet states, there is a 
corresponding probability associated with a temperature 
T. The effective coupling J depends on which state is 
chosen for the two spins. Energy, corresponding effective 
coupling, and the probability to be in the singlet/triplet 
states are listed in Table Hf. Note that if two spins are in 
Itt) or 14-4-) state, the effective coupling they mediate has 
opposite sign to the other cases. In the end, by using the 
RSRG-X process, a typical excited eigenstate associated 
with a temperature is approximated as the product state 
of these singlet and triplet pairs. One example of such a 
state generated by RSRG-X method for T 7 ^ 0 is depicted 
in Fig. [2] In this figure, singlet and triplet states are 
plotted in different colors. 

At each step of RSRG-X, effective coupling can be a 
negative or a positive number (J ss j-iTa depending on 
which of singlet or triplet states is chosen according to the 
Boltzmann distribution). After repeating RSRG-X steps 
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FIG. 1: A schematic plot of one step of real space renormal¬ 
ization group (RSRG/RSRG-X) procedure. Each solid line 
with coupling J connects two spins which are depicted by 
filled circles and the dashed line is an effective coupling J 
generated by RSRG/RSRG-X. 



FIG. 2: (Golor online) An example of spin configuration gen¬ 
erated by RSRG-X method to find approximate eigenstate 
of the XX spin chain corresponding to a non-zero tempera¬ 
ture. Every two spins are connected by a colored line: singlet 
(blue), triplet.|- 4 , (green), triplet.j..|. (cyan), triplet 4 , 4 , (red). 


many times, probability of having (positive or negative) 
smaller J’s increases and finally we have the following 
power law distribution function (we set the cutoff energy 
to be 1 ): 


P{J) 


f|J|-i+“ ,if|J|<l 
0 , otherwise 


( 11 ) 


with a approaching zero slowly with decreasing energy 
scale, or disorder strength increases with lowering energy 
scale. This is the same behavior as in RSRG, which es¬ 
tablishes their asymptotically exact nature. We choose 
this distribution function in our exact numerical calcula¬ 
tions. In this distribution function small a corresponds to 
strong disorder regime and large a corresponds to weak 
disorder regime. 

In the fermionic representation of the Hamiltonian, Eq. 
we need to know the distribution of fermions in en¬ 
ergy levels corresponding to the eigenstate obtained in 
RSRG/RSRG-X method. Eirst we calculate the single 
particle energy levels by diagonalizing the Hamiltonian 
in Eq. ([S]). Due to the particle-hole symmetry of the 
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Hamiltonian, energy levels form ±_B pairs. For a system 
with N sites, there are N/2 pairs of energy levels (assum¬ 
ing even N). Each pair of energy levels corresponds to an 
effective bond in the RSRG/RSRG-X method. For the 
ground state of the Hamiltonian, for each step of RSRG, 
we put the two spins coupled with the maximum J to 
be in the ground state and correspondingly we put one 
fermion in the negative sector of each pair of ±E and thus 
all negative energy levels are filled. The procedure is dif¬ 
ferent in the RSRG-X application for excited eigenstates. 
At each step of RSRG-X, number of fermions (0, 1, or 2) 
and which energy level(s) we put them in, depend on the 
state of the two spins (the two spins can be in singlet or 
either of triplet). For N/2 pairs of energy levels, we have 
N/2 RSRG-X steps; for nth step, the corresponding pair 
of energy levels are the nth and (N — n + l)th energy 
levels. We put no fermion if two spins are chosen to be 
in 14-4-) ■ If Boltzmann distribution chooses two spins to 
be in |tt)) then we put two fermions, one in each of two 
specified energy levels. For the case of ■^[|t4') i lit)]; we 
put one fermion in nth or (TV — n + l)th level depending 
the randomly chosen state is a ground state or an excited 
state (by considering the sign of the Jmax)- 


IV. RESULTS 

With regard to entanglement, predictions of the 
RSRG-X method are the following. The excited eigen¬ 
state is approximated as the product state of the singlet 
and triplet states obtained in all steps. When the system 
is divided to two subsystems, it may cut singlet and/or 
triplet states. For each singlet and triplet^j, crossing the 
boundary, there is one entanglement energy e at exactly 
0 and the rest are ±oo. In addition, corresponding Klich 
eigenmode for each e = 0 is highly localized at two spins 
in singlet or triplet-fj,. But since spins in triplet-i-j- and 
triplet^j, state are not entangled and do not contribute 
to the entanglement, Klich eigenmodes do not display the 
location of spins in these two states. 

In this section we analyze the predictions of the RSRG- 
X method. In first subsection, by looking at the entan¬ 
gled modes, we explain how they represent the entan¬ 
gled spins. In the subsequent subsection, we examine 
the predictions of RSRG-X about entanglement energies. 
For both subsections, we distinguish between strong and 
weak disorder limits to see in which limit RSRG-X pre¬ 
dictions are accurate. 


A. Entangled modes 

In our previous papers^ we explained that MEM re¬ 
veals important physics about the system. For example 
MEM displays location of the pair of spins in singlet state 
that crosses the boundary of subsystem (and thus con¬ 
tribute to the entanglement) in the ground state of the 


XX chain. Here, we show that MEM displays the loca¬ 
tion of spins that are in singlet or triplet-f-j, that cross the 
boundary for the excited states of the XX chain. 

In Fig. 131 there are four subplots, by which we want 
to demonstrate how MEM displays the location of en¬ 
tangled spins. In all of them we set TV = 60 and T = 1. 
There are three panels for each subplot. First panel is 
the presentation of the RSRG-X where we show each two 
spins in the singlet/triplet states connected by a colored 
line, second and third panels are the MEM associated 
with the specified partition of the system shown in the 
plot by a vertical line. In subplot (a), first panel, we 
show the RSRG-X of a sample that has only one singlet 
crossing the boundary and we choose a = 0.1 (strong dis¬ 
order regime). In subplot (b), we choose a sample that 
has only one triplet-i-j, crossing the boundary and also 
with a = 0.1. In second panel of both plots we plot the 
corresponding MEM. We repeated the same calculations 
in subplots (c) and (d), but instead we choose a = 0.9 
(weak disorder regime). 

Gomparing two cases of strong and weak disorder, we 
observe that: First, for the strong disorder case MEM 
is highly localized at the location of spins that are pre¬ 
dicted by RSRG-X to form singlet or triplet-j-j., clearly 
demonstrating its accuracy. For the weak disorder case, 
while MEM is more extended, we still see that it is cen¬ 
tered around spins that RSRG-X predicts to form singlet 
or tripletij, that cross the boundary. In this case what 
get entangled are effective spins renormalized through 
RSRG-X procedures, and MEM gives their profiles. 

Second, to see that these effective spins are intrinsic of 
the system and independent of the boundary, we choose 
a different partition in the third panel of each of the sub¬ 
plots and we see that still MEM is localized at the same 
entangled (effective) spins. More quantitatively, we cal¬ 
culate overlap between MEM’s corresponding to different 
subsystem sizes in Table iHl and IHII We see that MEM’s 
of case with different subsystem sizes but correspond¬ 
ing to same entanglement configuration are almost iden¬ 
tical, although in the week disorder case their overlaps 
are slightly smaller than the strong disorder case. Thus, 
MEM is a property of the system and it is independent 
of subsystem boundary. 

It is possible to have cases in which multiple sin¬ 
glet/triplet states cross the boundary of subsystem. We 
present one example in the following. Eig. |4]is a case with 
one singlet, one triplet-j-j,, one triplet-i-i-, and one triplet 4 . 4 , 
cross the boundary, and thus two entangled pairs of spins 
cross the boundary. In this figure, in which there are two 
entangled modes, we provide two Klich eigenmodes in 
second and third panels, corresponding to the two small¬ 
est magnitude of e’s. Both of them are localized at the en¬ 
tangled spins (singlet and triplet't- 4 .) and we also see that 
they do not display the location of the spins in triplet 4 - 4 ' 
or triplet 4 . 4 . since these two states do not contribute to 
the entanglement. 

Also, we find that in the strong disorder regime where 
MEM is localized at two entangled points, sign of the 
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FIG. 3: (color online) Demonstrating how MEM displays the location of entangled spins. In all four plots we set = 60 and 
T = 1. Subplots (a) and (b) correspond to the strong disorder {a = 0.1); (c) and (d) correspond to the weak disorder (a = 0.9). 
In subplots (a) and (c) there is only one entangled singlet (blue bond) and in subplots (b) and (d) there is only one entangled 
triplet.]- 4 _ (green bond). First panel of each subplot displays the RSRG-X method. Each two spins in singlet/triplet states are 
connected by a line with a specific color. Singlet (blue), tripleti. 4 . (green), triplet-j-j- (cyan), and triplet 4 . 4 . (red). Second and 
third panels show the MEM with specified subsystem size Na- 


TABLE II: Overlap between MEM’s corresponding to differ¬ 
ent subsystem sizes for a — 0.1 and a = 0.9, same setting as 
Fig. (3] subplots (a) and (c) respectively. Subsystem size with 
A = 55 corresponds to no entanglement. Other subsystem 
sizes correspond to the same case of one singlet crossing the 
boundary. 


a = 

- 0.1 

10 

1 — 1 

II 

Na = 19 

CO 

0 

0 

II 

Na 

= 

55 

Na 

= 15 

1 

1 

1 

1 

1 

X 

10’ 

-ib 

Na 

= 19 

1 

1 

1 

1 

1 

X 

10' 

-16 

Na 

= 30 

1 

1 

1 

1 

1 

X 

10 

-16 

Na 

= 40 

1 

1 

1 

1 

1 

X 

10' 

-16 

Na 

= 55 

1 X 10“i® 

1 X 10"^® 

1 X 10”1® 

1 X 10"i® 



1 



a = 0.9 Na = 1^ Na = 19 Na = 30 Na = 40 Na = 55 


Na 

= 15 

1 

0.99 

0.99 

0.99 

0.0004 

Na 

= 19 

0.99 

1 

0.99 

0.99 

0.0004 

Na 

= 30 

0.99 

0.99 

1 

0.99 

0.0004 

Na 

= 40 

0.99 

0.99 

0.99 

1 

0.0004 

Na 

= 55 

0.0004 

0.0004 

0.0004 

0.0004 

1 


MEM at these two points have a physical information. 
By looking at different examples we find the following 
rule. Remembering that each singlet/triplet contains 
specific number of fermions by Eq. m, we count num¬ 
bers of fermions by all singlet and triplet included inside 


TABLE III: Overlap between MEM’s corresponding to differ¬ 
ent subsystem sizes for a = 0.1 and a = 0.9, same setting 
as Fig. (3] subplots (b) and (d) respectively. Subsystem size 
with N = 5 corresponds to no entanglement. Other subsys¬ 
tem sizes correspond to the same case of one triplet 4 - 4 . crossing 
the boundary. 


a = C 

1.1 Na = 5 

to 

0 

Na = 30 

Na = 33 

1- 

CO 

II 

Na = 

■■ 5 

1 

1 X 10"^^ 

1 X 10"“ 

1 X 10"“ 

1 X 10"“ 

Na = 

20 1 

X 10"“ 

1 

0.99 

0.99 

0.99 

Na = 

30 1 

X 10““ 

0.99 

1 

1 

0.99 

Na = 

33 1 

X 10”“ 

0.99 

1 

1 

0.99 

Na = 

55 1 

X 10““ 

0.99 

0.99 

0.99 

1 

a = 

= 0.9 

Na=5 

Na = 20 TVa = 30 Na = 33 Na ^ 37 

Na 

= 5 

1 

0.004 

0.002 

0.002 

0.002 

Na 

= 20 

0.004 

1 

0.86 

0.83 

0.83 

Na 

= 30 

0.002 

0.86 

1 

0.83 

0.83 

Na 

= 33 

0.002 

0.83 

0.83 

1 

0.97 

Na 

= 55 

0.002 

0.83 

0.83 

0.97 

1 


the entangled singlet or triplet 4 - 4 , . In the case of en¬ 
tangled singlet if number of fermions is odd, then both 
non-zero magnitude of MEM have same sign, otherwise 
they have opposite signs. The rule is opposite for the case 
of one entangled triplet 4 - 4 ,. For example, in Fig. HI inside 
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FIG. 4: (color online) An example of having one singlet (blue), 
one tripletj^ 4 . (green), one triplet-|-|- (cyan), and one triplet 4 , 4 . 
(red) cross the boundary. First panel is the configuration of 
spins making singlet/triplet generated by RSRG-X method. 
Second and third panels are the Klich eigenmode correspond¬ 
ing to the two smallest entanglement energy e'^. 


the entangled singlet there are 0 singlet, 4 triplet 4 ' 4 ,s, 1 
triplet-i-i-, and 4 triplet 4 , 4 , 8 , so number of fermions included 
in is 6 and signs of the Klich eigenmode at the location of 
two entangled singlet are opposite. Inside the entangled 
triplet-|- 4 , there are 2 singlets, 7 triplet'|' 4 ,s, 5 triplet-i-i-s, and 
6 triplet 4 , 4 , 8 , so number of fermions included in is 19 and 
signs of the Klich eigenmode at the location of two en¬ 
tangled triplet-|- 4 , are opposite, which are consistent with 
the proposed rule. 

To quantify how much MEM is spread over sites, we 
use inverse participation ratio (IPR): 


IPR = 


1 

ETW 


( 12 ) 


If IPR is close to 2, it means that MEM is very local¬ 
ized at the positions of two spins, and when IPR is big¬ 
ger than 2, there is an spread of the MEM over a few 
sites and IPR/2 is the spatial extent of effective spins. 
In Pig. [5l we plot IPR of the MEM as a changes for 
two different cases of having only on singlet or only one 
triplet-j- 4 , crossing boundary for a non zero temperature 
T = 0.1. For both of them we see that in the strong dis¬ 
order regime (small a) IPR is close to 2 and this means 
that MEM is very localized at the position of two entan¬ 
gled spins. This is consistent with RSRG-X predictions, 
since we know that in this method (which is asymptot¬ 
ically exact in the strong disorder regime) at each step 
two spins make a singlet or triplet. As we approach to 
the weak disorder regime, IPR increases; meaning that 
MEM is spread over a few spins and effective spins form 
entanglement. 

Also, to show that effective spin size is independent 
from the system size, we provide a plot of IPR of the 
MEM versus system size iV (we set Na = N12), for two 
cases of strong and weak disorder regime when we have 
one singlet or one triplet'i- 4 , crossing the boundary in Fig 
. El Beside small fluctuations, we see that IPR is con¬ 
stant as system size N changes. This is true for having 
one singlet or one triplet'i- 4 , crossing the boundary, i.e the 



FIG. 5; (color online) IPR of the MEM as a changes for a 
system with N = 100, Aa = 50, T = 0.1. Blue line corre¬ 
sponds to the case of having only one entangled singlet and 
the red line corresponds to one entangled triplet 4 . 4 ,. IPR be¬ 
comes bigger as a increases for both cases. For each a, IPR 
is calculated for 40 different random J’s and for each of J’s 
distributions, thermal averaging is done over 20 excited states 
generated by RSRG-X corresponding to the specific tempera¬ 
ture T = 0.1. Overall, there are 800 number of data for each 
a. 



FIG. 6 : (color online) IPR of the MEM as N changes between 
100 and 300 for two cases of a low temperature (T = 0.01) and 
a high temperature (T = 1). Blue line corresponds to a = 
0.1 and red line corresponds to a = 0.9. Left (right) panels 
correspond to one entangled singlet (one entangled triplet 4 - 4 ,) 
crossing the boundary. Despite small fluctuations, IPR stays 
constant as N changes, although for the case of weak disorder, 
fluctuations are bigger. For each N, IPR is calculated for 40 
different random J’s and for each of J’s distributions, thermal 
averaging is done over 20 excited states generated by RSRG- 
X corresponding to the specific temperature. Overall, there 
are 800 number of data for each N. 


effective spin is a property of the of the MEM. 

More importantly, we consider the change of IPR of the 
MEM as temperature T changes to see how approach¬ 
ing excited states affects the validity of the RSRG-X 
method. In Fig. [7] we plot IPR of the MEM for a range 
of temperatures, from very small T = 0.001 to a large 
T = 100, again for two cases of having only one singlet 
or one triplet-j- 4 , crossing the boundary. We see that in 
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FIG. 7: (color online) IPR of the MEM as temperature 
changes between 0.001 and 100 in log-linear scale for two 
case of a = 0.1 (blue line) and a = 0.9 (red line). We set 
N = 100, Na = 50 for all calculations. In left panel, cases 
with only one entangled singlet and in right panel cases with 
only one triplet-t- 4 . are considered. For each T, IPR is cal¬ 
culated for 40 different random J’s and for each of J’s dis¬ 
tributions, thermal averaging is done over 20 excited states 
generated by RSRG-X corresponding to the specihc temper¬ 
ature. Overall, there are 800 number of data for each T. 


both cases IPR is close to 2 for strong disorder case, i.e. 
MEM is highly localized at the locations of two entan¬ 
gled spins even for a highly excited state associated with 
a high T; this is consistent with the predictions of the 
RSRG-X method for a finite temperature and thus in 
the strong disorder regime RSRG-X method is still valid 
for excited states corresponding to T 7 ^ 0 (for the ground 
state, T = 0, we concluded to the same conclusion in our 
previous paper—). On the other hand, for the weak dis¬ 
order case, IPR is always bigger than 2 and it becomes 
bigger for higher temperatures, i.e. spatial extent of the 
effective spins increases as temperature increases. 


Fig. | 8 ] illustrates the behavior of the MEM with vary¬ 
ing bond length. In this figure, IPR of the MEM of 
samples with only one singlet or one triplet-i-j. crossing 
the boundary are considered. IPR of the MEM for two 
temperatures T = 0.01,1 and three disorder strengths 
a = 0.6, 0.8, 0.9 is calculated. We see that IPR initially 
increases, followed by a decrease, and then saturates to 
a constant for long bond lengths, with the constant de¬ 
pending on disorder strength. This is consistent with the 
expectation that the system flows to a fixed point under 
RSRG-X, as initially it is the bare spins that form singlet 
or triplet pairs (hence the small IPR) while renormaliza¬ 
tion toward effective spins (with growing sizes) occur in 
the meantime; eventually the renomalization process sat¬ 
urates, and the size of the effective spins (measured by 
IPR of MEM) no longer depends on the bond length. 
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FIG. 8 : (color online) IPR of the MEM versus bond length 
when system size goes between 100 and 1500 with step of 
100. IPR of the MEM of samples with only one singlet or 
one triplet-t' 4 . crossing the boundary are considered. IPR is 
calculated for two temperatures T = 0.01 (left panel) and 
T = 1 (right panel) and in each panel for three different 
a = 0.6, 0.8, 0.9. For each N, IPR is calculated for 40 different 
random J’s corresponding to a specific a and for each of J’s 
distributions, thermal averaging is done over 20 excited states 
generated by RSRG-X corresponding to the specific temper¬ 
ature. Range of bond length which is between 1 and 1500 is 
divided to equal segments and for each segment average of the 
IPR is plotted. For both small and large temperatures IPR 
approaches to a constant dependent on the disorder strength 
a. 


B. Entanglement Energy 

In the RSRG-X treatment, two spins coupled by the 
strongest bond at a particular step of RG form a sin¬ 
glet or one of the three triplets. For each entangled sin¬ 
glet or triplet-j-j,, corresponding entanglement energy is 
exactly zero. Otherwise entanglement energies are ±00 
(i.e., no entanglement). In the following, first we exam¬ 
ine this prediction of RSRG-X by looking at the entan¬ 
glement energies of a specific excited state of one sam¬ 
ple (without doing thermal or sample averaging) with 
N = 100, Na = 50 at T = 1, that has two pairs of entan¬ 
gled spins: on singlet and one triplet-i-j,. Table ITVl and IVl 
show the entanglement energies e^’s and the correspond¬ 
ing for strong and weak disorder cases respectively. 
In the strong disorder case {a = O.I), two e’s are very 
close to zero and beside a few finite e’s, many others are 
±c». However, for weaker disorder case {a = 0.9), two 
smallest e’s are distinguishable from zero and also there 
are more finite e’s. We see that in the strong disorder 
case, entanglement energies are in good agreement with 
the predictions of RSRG-X. 

To see the overall behaviour of smallest entanglement 
energy as disorder strength changes, we plot average of 
the smallest e magnitude versus a in Fig. [5] We set 
N = 100, = 50, T = 0.1 and select samples with 

one entangled singlet or one entangled triplet^j,. We see 
that in the strong disorder regime, smallest e is close to 
zero and increases as a increases. This is consistent with 
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FIG. 9: (color online) Average of the smallest e magnitude 
versus disorder strength. We set N — 100, Aa = 50, T = 0.1 
and select samples with one entangled singlet or one entan¬ 
gled triplet.|. 4 .. For each a, e is calculated for 80 different 
random J’s and for each of J’s distributions, thermal aver¬ 
aging is done over 40 excited states generated by RSRG-X 
corresponding to the specific temperature T = 0.1. Overall, 
there are 3200 number of data for each a. 

the expectation that accuracy of RSRG-X decreases with 
decreasing randomness strength, just like RSRG. 

To see how the accuracy depends on tempera¬ 
ture/energy density, we plot the average of the smallest e 
magnitude versus temperature T for two cases of strong 
and weak disorder in Fig. [TUI We set N = 100, Na = 50 
and consider samples with one entangled singlet or one 
entangled tripletfi. We see that in the strong disorder 
case {a = 0.1), e increases very slowly with temperature, 
and remains very small for the entire temperature, attest¬ 
ing to the accuracy of RSRG-X. For the weaker disorder 
case (a = 0.9), e is close to zero for low temperatures and 
increases considerably faster with increasing T. This in¬ 
dicates the quantitative error of RSRG-X increases with 
temperature/energy density. 

V. CONCLUDING REMARKS 

We have studied in this paper entanglement energies 
and entangled modes of highly excited states in a ran¬ 
dom one-dimensional XX chain, by mapping the system 
onto a free fermion model. We found qualitative agree¬ 
ment with predictions based on the newly developed real 
space renormalization group method for excited states 
(RSRG-X)ii, lending support to its validity. Further¬ 
more we found RSRG-X is quantitatively accurate in 


the strong disorder regime. As in the case of RSRG for 
ground state there are quantitative errors for weaker ran¬ 
domness, and such error grows with increasing tempera¬ 
ture/excitation energy density. This is not surprising as 
there are resonances between nearly degenerate excited 
states not included in the RSRG-X approximation; this 
source of error does not exist for ground state. Our work 
is complementary to an earlier studyi, which focused on 



FIG. 10: (color online) Average of the smallest e magnitude 
versus temperature T. We set N = 100, Na = 50 and consider 
samples with one entangled singlet or one entangled triplet^-j.. 
For each T, e is calculated for 80 different random J’s and 
for each of J’s distributions, thermal averaging is done over 
40 excited states generated by RSRG-X corresponding to the 
specific temperature. Overall, there are 3200 number of data 
for each T. 

entanglement entropy of the same model. As we demon¬ 
strated here and in earlier works^i^, entanglement ener¬ 
gies and in particular, entangled modes provide valuable 
additional information about the system. 

In a broader context, our work is part of on-going at¬ 
tempt to study entanglement properties of highly excited 
states. Such properties have been attracting consider¬ 
able attention recently, especially in the context of ther- 
malization and many-body localization^^. While most of 
the existing studies focus on the scaling behavior entan¬ 
glement entropy, our work suggests other entanglement 
measures also deserve investigation. 
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